Theoretical study of the experimentally important states of dysprosium 
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Configuration interaction method is used to calculate transition amplitudes and other properties 
of the low states of dysprosium which are used in cooling and in study of the time variation of 
the fine structure constant and violation of fundamental symmetries. The branching ratio for the 
cooling state to decay to states other than ground state is found to be smaller than 10~*. The 
matrix element of the weak interaction between degenerate states &t E — 19797.96 cm~^ is about 
2 Hz which is consistent with the experimental limit \Hw\ ~ |2.3± 2.9(statistical) ±0.7(systematic)| 
Hz [A. T. Nguyen, D. Budker, D. DeMille, and M. Zolotorev, Phys. Rev. A 56, 3453 (1997)] 
and points to feasibility of its experimental measurement. Applications include search for physics 
beyond the standard model using the PNC isotopic chain approach. 

PACS numbers: 31. 15. am, 32.70.Cs, 31.30.jg 
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I. INTRODUCTION 

Atomic dysprosium is proved to be a valuable object 
to study fundamental problems of modern physics. It 
has been recently used in search for the time variation 
of the fine structure constant [l| , study of possible vari- 
ation of the fine structure constant due to variation of 
the gravitation field |2| , measurement of the parity non- 
conservation (PNC) [3], etc. Recent progress in trapping 
and cooling of dysprosium atoms [^-l6| opens new excit- 
ing possibilities. 

The most interesting feature of dysprosium which has 
inspired its use in the study of PNC and variation of 
the fine structure constant is the existence of two almost 
degenerate states of the same total momentum and op- 
posite parity. Energy interval is so small that its actual 
value is determined by hyperfine structure and isotope 
shift. This leads to strong enhancement of both effects. 
The states are at the energy E — 19797.96 cm~^ and 
both have total momentum J — 10. Following Ref. 
we use notation A for the even state and notation B for 
the odd state. 

Another state of great interest is the odd state e.t E — 
23736.60 cm~^ which is used in cooling. We use notation 
C for this state. Our present work mostly focuses on these 
three states while some other states are also considered. 

The use of the rare-earth atoms, including dysprosium, 
in atomic PNC study was first suggested in Ref. [8] , and 
for the search of the time-variation of the fundamental 
constants in Refs. [1, • Dysprosium was studied the- 
oretically in our previous works The work of 
Ref. [13, finks the change of frequency of the tran- 
sition between states A and B to the time-variation of 
the fine structure constant. It was used in Ref. for 
the interpretation of the measurements. In Ref. ll| we 
calculated the matrix element of the spin-independent 
parity-violating weak interaction between states A and 
B. 

Dysprosium has many stable isotopes and is a good 
candidate to study PNC ratio for isotope chains. Such 
study may reveal important information on the physics 



beyond the standard model Accurate atomic cal- 

culations of the PNC effect are not needed for this 
study. However, reliable estimation of the effect is im- 
portant to determine the feasibility of the measurements. 
The result of our previous calculations is (A|iy|B) = 
70(40) Hz [ll|. Later measurements [H lead to the limit 
(A|iy|B) = |2.3 ± 2. 9 (statistical) ± 0.7(systematic). Al- 
though the measured value is not in strong disagreement 
with theoretical prediction, given the large uncertainty 
of the latter, the experiment did not confirm the large 
PNC effect that was hoped for in this system. 

In present paper we revisited the PNC calculations and 
found that inclusion of more configurations pull the value 
of the weak matrix element down to about 2 Hz which 
is consistent with the measurements. The small value 
of the matrix element is the result of strong cancella- 
tion of different contributions. Dominant contributions 
are larger than the final result by more than order of 
magnitude. This means that further cancellation to even 
smaller number is highly unlikely and the measurements 
might be possible on about the same level of sensitivity 
which has been already achieved in Ref. 

Another motivation for this work is due to dyspro- 
sium cooling at Berkeley [1, Q and Urbana [H]- We 
study the cooling state C to find the branching ratio of 
the transitions from this state to the states other than 
the ground state and to the ground state. High value 
(;» 10^"*) for this ratio would be a problem for cooling. 
Our present calculations show that the branching ratio 
is in fact smaller than lO^'*. 



II. METHOD 

In present work we use the version of the configuration 
interaction (CI) method which was first developed for 
iron atom [15[ and then used for other many-electron 
atoms including dysprosium [ll] . See these works for the 
detailed discussion. 

The effective Hamiltonian for Ny valence electrons 
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TABLE L Configurations and effective core polarizabilities ap 
(a.u.) used in tfie calculations. 



N 


Parity 


Configuration 


Op 


1 


Even 


4/^"6s^ 


0.4 


2 


Even 


4/^°6s5d 


0.4006 


3 


Even 


4/^6s^6p 


0.4039 


4 


Even 


4/^5d6s6p 


0.389 


5 


Even 


4/ "6p 


0.4 


6 


Even 




0.4 


7 


Odd 




0.3947 


8 


Odd 




0.3994 


9 


Odd 


4:f°6s6p 


0.397 


10 


Odd 




0.4 


11 


Odd 


4f5d6p'^ 


0.4 


12 


Odd 


'if6s6p'^ 


0.4 



(— Q;p/2r^). We treat ap as fitting parameters. The val- 
ues of ttp for each configuration of interest are presented 
in Table HI They are chosen to fit the experimental posi- 
tion of the configurations relative to each other. For all 
configurations the values of ap are very close. This is not 
a surprise since the core is the same for every configura- 
tion. Small difference in ap for different configurations 
simulates the effect of incompleteness of the basis and 
other imperfections in the calculations. 

To calculate electric dipole transition amplitudes, life- 
times, magnetic dipole and electric quadrupole hyperfine 
structure constants we use the time-dependent Hartree- 
Fock method (equivalent of the random phase approxi- 
mation) combined with the CI technique, see Ref. [3| 
for a detailed discussion. 



(A'^„=12 for dysprosium) has the form 

iJ^ff = ^/ii,-^^e>,„ (1) 

i=i i<j 

hi{ri) is the one-electron part of the Hamiltonian 

hi = ca--p^ [13 - l)mc^ -I- Vcore + SV. (2) 

r 

Here a and j3 are Dirac matrixes, Vcore is Hartree-Fock 
potential due to core electrons and 5V is the term which 
simulates the effect of the correlations between core and 
valence electrons. It is often called polarization potential 
and has the form 



Here ap is polarization of the core and a is a cut-off 
parameter (we use a — ub)- 

Table H] lists configurations considered in present work. 
The self-consistent Hartree-Fock procedure is done for ev- 
ery configuration separately. Then valence states found 
in the Hartree-Fock calculations are used as basis states 
for the CI calculations. It is important for the CI method 
that the atomic core remains the same for all configura- 
tions. We use the core which corresponds to the ground 
state configuration. Change in the core due to change 
of the valence state is small and can be neglected. This 
is because core states are not sensitive to the potential 
from the electrons which are on large distances (like 6s, 
6p and 5d electrons). The 4/ electrons are on smaller dis- 
tances and have larger effect on atomic core. However, 
in all the cases (see Table |l| only one among about ten 
4/ electrons change its state. Therefore their effect on 
atomic core is also small. More detailed discussion on the 
effect of valence electrons on atomic core can be found in 
Refs. 

The form of the SV in ([3]) is chosen to coincide with 
the standard polarization potential on large distances 



III. RESULTS AND DISCUSSION 

Table HIl presents the results of calculations for the en- 
ergies, g-factors, magnetic dipole and electric quadrupole 
hyperfine structure constants and lifetimes of some low 
states of dysprosium. This includes the degenerate states 
A and B, cooling state C and other odd states D to G 
which may also present an interest for cooling and quan- 
tum information processing [l^. Calculated values are 
compared with available experimental data. Good agree- 
ment for the energies is mostly due to the fitting. Theo- 
retical uncertainty for the hyperfine structure constants 
is on the level of 20 to 30%. The uncertainty for lifetimes 
is determined by uncertainties for electric dipole transi- 
tion amplitudes. The uncertainties for the amplitudes 
is also on the level of 20 to 30% with the exception of 
the extremely small amplitudes (<C 1 a.u.) where uncer- 
tainty might be higher. The amplitudes are presented in 
Table IIIII Maximum disagreement between theory and 
experiment in Table [H] is for the lifetime of the state A 
where calculated value is two times larger than the exper- 
imental one. The lifetimes of the most of the states are 
sensitive to the mixing of the 4/^°6s6p with other odd 
configurations. For example, the experimental lifetime of 
the state C is reproduced in the calculations when state 
C is a pure state of the 4/^°6s6p configuration with very 
little admixture of other configurations. This has been 
achieved by properly choosing the fitting parameters ap 
(see Table HI). 

Table Hill shows decay channels for the states listed in 
Table [Hi The data is based on the calculated electric 
dipole transition amplitudes and experimental energies. 
No electric quadrupole or magnetic dipole amplitudes 
were taken into account. Test calculations show that 
their contribution is negligible. The analysis of the data 
for state A shows that it decays mostly to two states of 
the 4/^^656^ configuration. This is because the leading 
configuration of state A is the 4/^°6s5(i configuration and 
an electric dipole transition between the two states can 
be reduced to the allowed single-electron 6p 5d transi- 
tion. In contrast, transitions to the states of the 4/^5(i6s^ 
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TABLE II: Experimental and theoretical energies, g-factors, hyperfine structure constants and lifetimes for some low-lying 
states of ^^^Dy. For theoretical uncertainties see discussion in the text. 



N 


leading 


Term 


J 


Energies ( 


cm'i) 


3-factors A (MHz) 


B (MHz) 


Lifetime 




config. 






Expt." 


Calc. 


Expt." 


Calc. Expt. 


Calc. 


Expt. Calc. 


Expt. 


Calc. 












States of special interest 










GS' 




^I 


8 


0.00 





1.24 


1.24 163" 


160 


1153" 1193 






A 


4/i°5d6s 


'[10] 


10 


19798 


19786 


1.21 


1.21 159^* 


140 


1865'' 1629 


7.9 ^is'' 


16 fis 


B 


4/^5d^6s 


10 


19798 


19770 


1.367 


1.368 218'' 


202 


2060'' 2413 


> 200 fis"^ 


0.14 s 


C 


4/^°6s6p 


(8,ir 


9 


23737 


25200 


1.22 


1.22 122" 


136 


1842" 2096 


4.8 ns'' 


4.7 ns 














Some other states 










D 


4/^5d6s2 




8 


7566 


7563 


1.35 


1.35 


131 


1901 




6.9 ms 


E 


4/^5^65^ 


7jo 


9 


9991 


9944 


1.32 


1.32 


125 


2901 




3 ms 


F 






9 


13496 


14634 


1.23 


1.23 


144 


4000 




21 ^s 


G 


4/^°6s6p 


(8,2)° 


9 


17727 


18092 


1.25 


1.26 


194 


1582 


2 fis'' 


2.9 ^is 



"Reference i22i 
''Ground state 
'^Reference |23l425l | 
Reference 01 
'^Reference [21| 



TABLE HI: Decay channels for states A - G from Table HIl 



Upper 




Lower state 




Energy 




K/I|D||*>I 


Probability 


state 


Config. 


Term 


J 


cm^^ 


cm^ 


(a.u.) 


s-i 


A 


4/^5d6s^ 


'K° 


10 


12893 


6905 


0.056 


0.982E+02 




4/i°6s6p 


(8,2)'^ 


10 


17513 


2285 


0.276 


0.874E+02 




4f5d6s'^ 


7p 


9 


9991 


9807 


0.049 


0.218E+03 




4f5d6s'^ 




9 


13496 


6302 


0.192 


0.887E+03 




4/i°6s6p 


(8,1)° 


9 


15972 


3826 


2.34 


0.296E+05 




4f5d6s'^ 


7j^o 


9 


16717 


3080 


0.441 


0.547E+03 




4/i°6s6p 


(8,2)° 


9 


17727 


2071 


5.86 


0.294E+05 


B 


4/i°5d6s 


'[9] 


10 


18463 


1335 


0.069 


0.109E+01 




4/i°5d6s 


'[8] 


9 


17515 


2283 


0.073 


0.608E+01 




4/i°5d6s 


'[9] 


9 


19241 


557 


0.0045 


0.339E-03 


C 


4/i06s^ 




8 





23737 


12.28 


0.215E+09 




4/i°5d6s 


'[7] 


8 


17613 


6123 


0.221 


0.120E+04 




4/i°5d6s 


'[8] 


8 


18903 


4833 


0.265 


0.845E+03 




4f%s^ 


^K2 


8 


19019 


4717 


0.057 


0.362E+02 




4f°5d6s 


'[9] 


8 


20194 


3543 


0.036 


0.601E+01 




4f6s6p 


CV2,V2) 


8 


20790 


2947 


0.099 


0.265E+02 




4/i°5d6s 


'[7] 


8 


21603 


2134 


0.226 


0.530E+02 




4f°5d6s 


'[8] 


9 


17515 


6222 


0.411 


0.434E+04 




4/i°5d6s 


'[9] 


9 


19241 


4496 


0.122 


0.144E+03 




4/i°5d6s 


'[10] 
M9] 


9 


20209 


3528 


0.222 


0.231E+03 




4/i°5d6s 


9 


22046 


1691 


0.659 


0.224E+03 




4/i°5d6s 




9 


23218 


518 


0.0076 


0.846E+03 




4/i°5d6s 


'[9] 


10 


18462 


5274 


0.502 


0.394E+04 




4/i°5d6s 


'[10] 


10 


19798 


3939 


0.467 


0.142E+04 




4/^°5d6s 


Mio] 


10 


22487 


1249 


1.039 


0.225E+03 


D 


4/i06s2 




8 





7566 


0.053 


0.146E+03 




4f%s^ 




7 


4134 


3431 


0.0017 


0.136E-01 


E 


4f%s^ 




8 





9991 


0.059 


0.369E+03 


F 


4f%s^ 




8 





13496 


0.424 


0.471E+05 


G 


4f%s^ 


'"1 


8 





17727 


0.897 


0.478E+06 
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configuration can only go via configuration mixing. The 
calculated lifetime of the state A is larger than the ex- 
periment. This means that the transition amplitudes are 
smaller. It is unlikely that the amplitudes of the tran- 
sitions between state A and states of the 4/^*^656^ con- 
figuration are underestimated. This is because the latter 
states are almost pure, with only about 3% admixture of 
other configurations. It is more likely that some of the 
other amplitudes are underestimated because of too small 
mixture of the Af^bdQs^ and 4/^°6s6j> configurations. 

The analysis of the decay channels of the cooling level 
C (see Table llll|) also indicates the sensitivity of the data 
to the mixing of the 4/^5(i6s^ and 4/^°6s6p configura- 
tions. Good agreement with experiment for the lifetime 
of the state C can only be achieved if this state is a pure 
state of the 4/^°6s6p configuration. On the other hand, 
this state is very close to the states of the 4/^5d6s^ con- 
figuration. Manipulating with the fitting parameters ap 
(Table H]) can easily lead to a situation when the states 
of the two configurations are strongly mixed and the life- 
time of the state C is larger than the experiment. There- 
fore, the lifetime of this state should be monitored in the 
fitting process. 

The data in Table Hill allows to estimate the ratio of the 
following probabilities: (a) decay of state C to all lower 
state other than the ground state, and (b) decay of the 
state C into the ground state. This ratio is (6±3) x 10~^. 
The fraction of atoms lost into metastable states must 
be even smaller because the number above includes all 
channels of decay without considering which of them end 
up in a metastable state. However, some lower states 
quickly decay to the ground state and do not produce 
any losses. A more detailed analysis would require a lot of 
extra work, however it is not needed since the branching 
ratio < 10""* is sufficiently good for the cooling Q. 

Table IIIII also shows the data for states D-G which 
might be useful for cooling and information process- 

ing ai- 

Our calculated value for the electric dipole reduced ma- 
tric element between states A and B is 0.024 ±0.010 a.u. 
This agrees well with the experimental value of 0.015(1) 
a.u. from Ref. [^. Note that this amplitude is zero in 
the non-relativistic limit. This is because states A and B 
have different spin (see Table |lTl and electric dipole op- 
erator cannot change it. In relativistic calculations the 
amplitude is not zero. However, it is small and this small 
value is the result of strong cancellations between differ- 
ent contributions. Table ITVl shows the largest in absolute 
values contributions to the electric dipole transition am- 
plitude between states A and B. All these contributions 
are due to the 4/ - 5d electric dipole transitions within 
the main configurations of states A and B. The sum is 
only 0.0017 a.u. which is more than an order of magni- 
tude smaller than the final answer. The final amplitude 
is the sum of many smaller contributions which contain 
all possible single-electron transitions. 

The situation is even more complicated for the ma- 
trix element of the weak interaction between states A 



TABLE IV: Contributions to the electric dipole transition am- 
plitude between states A and B (a.u.) 



Configurations 
State A State B 


Single-electron 
matrix element 


Partial 
contribution 


Sum 
(a.u.) 


4/""5d6s 4/"5d^6s 
4/"'5d6s 4f5d-^6s 
4f°5d6s 4f5d'^6s 


(4/5/2 |D|5d3/2> 
(4/5/2 |D|5d5/2> 
(4/7/2 |/3|5d5/2> 


-0.2634 
-0.0106 
0.2723 


-0.2634 
-0.2740 
-0.0017 



and B. Table |V] shows dominant contributions to this 
matrix element. There are strong cancellations between 
different contributions, mostly between terms contain- 
ing the (4/5/2 1 77,4' 1 5^5/2) and {6si/2\Hw\Qpi/2) single- 
electron matrix elements. The former of these matrix 
elements is small. It is not zero due to the effect of core 
polarization by the weak interaction: weak interaction 
acts on s and p core electrons changing the self-consistent 
Hartree-Fock potential which in turn leads to the 4/ - 5d 
transition between valence states. The (6s 1/2 6^1/2) 
integral is not small but its contribution is suppressed by 
small admixture of the configurations containing the 6s 
and 6p states to the main configurations. This makes the 
result very sensitive to configuration mixing. We found 
that inclusion of all configurations listed in Table U is im- 
portant for the weak matrix element. The effect of some 
configurations is indirect, via changing the coefficients of 
configuration mixing. For example, the sum of all largest 
contributions listed in Table IVl -2.44 Hz is very close to 
the final answer -2.13 Hz. This may make an impres- 
sion that the configurations not listed in Table |V] do not 
contribute. However, if e.g. the 4/^6s^6p configuration is 
not included, the weak matrix element becomes larger by 
more than four times. This configuration was missed in 
our previous calculations [Tl| which is probably one of the 
reasons of disagreement between theory and experiment. 
The result of Ref. corresponds to the situation when 
the weak matrix element between states A and B is domi- 
nated by the contribution of the (4/5/2 5^5/2) single- 
electron matrix element (first line of Table Ivp . This is 
in spite of the fact that most of the other configura- 
tions were also included. Probably incomplete inclusion 
of other configurations leads to underestimation of the 
configuration mixing which suppresses the contributions 
containing other single-electron matrix elements. 

Our present result for the weak matrix element 
(2 Hz) is consistent with the experimental value of 
\{A\Hw\B)\ = |2.3±2.9(statistical)±0.7(systematic)| d. 
This is about 20 times smaller than the values of indi- 
vidual contributions (see Table W} . Since there is no law 
of physics to make this matrix element exactly zero we 
believe that any further cancellation which would make 
the result even smaller than 2 Hz is highly unlikely. This 
probably means that the measurements of the PNG ef- 
fect can be done on about the same level of sensitivity 
which has been already achieved in Ref. 
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TABLE V: Contributions to the matrix element of the weak 
interaction between states A and B (Hz) 



Config 


urations 


Single-electron 


Partial 


Sum 


Left 


Right 


matrix element 


contribution 


(Hz) 



4/^°5d6s 4/^"5d6p 

4/i°5d6s Af°6s6p 

Af5d6s6p 4f5d6s'^ 

Af5d6s6p Af5d^6s 

4/^5d^6p 4/''5d^6s 



75/2 I 

{6si/2\Hw\6pi/2) -49.04 -10.25 

{5dy2\Hw\&P3/2} L60 -8.65 

{6pi/2\Hw\6s^/2) -1-35 -10.00 

{6ps/2\Hw\5d3/2} -4.28 -14.27 

{6pi/2\Hw\dsy2} 11.84 -2.44 
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